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ABSTRACT 

We present the discovery of a hot Jupiter transiting an F star in a close visual (0.3" sky projected 
angular separation) binary system. The dilution of the host star's light by the nearly equal magnitude 
stellar companion (~ 0.5 magnitudes fainter) significantly affects the derived planetary parameters, 
and if left uncorrected, leads to an underestimate of the radius and mass of the planet by 10% and 
60%, respectively. Other published exoplanets, which have not been observed with high-resolution 
imaging, could similarly have unresolved stellar companions and thus have incorrectly derived plane- 
tary parameters. Kepler-14b (KOI-98) has a period of P = 6.790 days and correcting for the dilution, 
has a mass of = 8.40lQi[gMj and a radius of Rp — 1.1361qq54 yielding a mean density of 
Pp = 7.1 ± l.lgcm"^. 

Subject headings: planetary systems — stars: individual (Kepler-14b, KIC 10264660, 
2MASS J19105011+4719589) techniques: spectroscopic, photometric 
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1. INTRODUCTION 

Kepler is a space-based telescope using transit pho- 
tometry to determine the frequ ency and characteris tics of 
planets and planetary systems ([Borucki et al.|[20lol ). The 
instrument was launched in March 2009 and is a wide 
field-of-view (115 square degrees) photometer comprised 
of a 0.95-meter effective aperture Schmidt telescope that 
monitors the brightness of about 150,000 stars. Recently, 
the first four months of photometric data were released 
and o ver 1200 trans i ting p lanet candidates were iden- 
tified (jBorucki et al.l 120111 ). Kepler-14b was identified 
among the 1235 candidates as KOT98. Because of its 
short period (P = 6.790 days) and its relatively deep 
transit signal, it was identified very early in the mission 
and quickly passed on to the Kepler Follow-up Program 
(KFOP) for further investigation. 

Kepler- 14b was scrutinized for evidence of astrophysi- 
cal false positives and survived the initial follow-up stage. 
We therefore gathered high resolution, high signal-to- 
noise ratio (SNR) spectra in order to extract precise ra- 
dial velocities (RVs) using the Fiber-fed Echelle Spectro- 
graph (FIES) on the Nordic Optical Telescope (NOT) at 
La Palma in October 2009. The observations yielded a 
spectroscopic orbit in phase with the photometric ob- 
servations from Kepler. Meanwhile, other ground based 
observations were gathered and speckle imaging revealed 
that the host star of Kepler- 14b was in fact not a single 
star, but a nearly equal magnitude binary system with a 
sky projected separation of only 0.29". This complicated 
the analysis of the system, and Kepler-14b was therefore 
not readied for publication together with the initial batch 
of five Kepler planets. 



2 



Buchhave et al. 



N 





-^nn ^ . I • . , , 'i ■ . , - I . , , I ■ . , ■ , I 

-4 -2 2 4 

Time since mid-lransil [hr] 



Figure 1. Light curve for Kepler-14b. The upper panel shows 
the photometry folded with the orbital period of 6.790 days. The 
fitted transit model is overplotted as a solid line. Residuals from 
the fit are displayed in the bottom panel. 



Table 1 

Speckle measurements of Kepler-14 



Band 


Separation" 


Position angle'' 


^mag 


V 


0.286 ±0.04 


143.67 ±0.07 


0.52 ±0.05 


R 


0.289 ±0.01 


143.67 ±0.07 


0.54 ±0.12 


I 


0.289 ± 0.02 


143.91 ±0.05 


0.45 ± 0.04 



In arcseconds. 
In degrees. 

In this paper, we present the confirmation of Kepler- 
14b as a planetary companion and analyze the photomet- 
ric and spectroscopic data taking into account the effect 
of the dilution by the stellar companion on the derived 
planetary parameters. Without the high spatial resolu- 
tion imaging, the stellar companion would not have been 
detected and Kepler-14 would have been published with 
incorrect planetary parameters. Other published transit- 
ing planets that have not been observed with high spatial 
resolution imaging might suffer the same problem. 

2. KEPLER PHOTOMETRY 
We use the long cadence photometry (29.4 minute ac- 
cumulations) of Kepler-14 (KIC 10264660, 19'^10'°50!12, 
+47°19'58'.'98, J2000, KIC r = 12.128 mag) obtained 
betw een 5 May 2009 to 7 September 2010 (QO through 
Q6) () Jenkins et al.ll2010bO. The photometrv is processed 
in an analysis pipeline penkins et al.ll2010c[ ) to produce 
corrected pixel data. Simple aperture photometry sums 
are formed to produce a photometri c time series, which 
was detrended as explained in Koc h et al.l (j2010t ). The 
photometric data folded with the orbital period of 6.790 
days are shown in Figure [TJ The numerical data are avail- 
able electronically from the Multi Mission Archive at the 
Space Telescope Science Institute (MAST) Web sitfE]- 



3. FOLLOW-UP OBSERVATIONS 

3.1. High-resolution speckle imaging 

Speckle observations of Kepler-14 were obtained on 6 
different nights between June and Oct 2010. The obser- 
vations were obtained with the dual channel WIYN tele- 
scope speckle camera recently described in iHorch et al.l 
(|2011f ). The data collection, reduction, and image recon- 
struction process is described in the aforementioned pa- 

^'^ http;/ /archive. stsci.edu/kepler/ 



Figure 2. Speckle image of Kepler-14 with its obvious companion 
star separated by 0.3"of nearly equal brightness (~ 0.5 magnitude 
fainter). 





Figure 3. J and Kg Palomar adaptive optics images of Kepler-14, 
showing a 2" X 2" field of view centered on the brighter star. 



per as well as in lHowell et al.l (|2011| ). The latter presents 
details of the 2010 season of observations for the Kepler 
follow-up program. 

A spatially close (0.3"), nearly equal brightness (^ 0.5 
magnitude fainter) companion star was easily noted in 
the reconstructed speckle images. Table [T] gives the 
weighted mean values for the separation, position angle, 
and magnitude difference for our six speckle observations. 
The observations are weighted by the native seeing dur- 
ing the time of the speckle data collection as determined 
by the data reduction routines when fitting known single 
point source speckle standard stars obtained near in time 
to the Kepler star observations. Figure [2] shows one of 
the reconstructed images of Kepler-14 with its obvious 
companion star. 

3.2. High-resolution Palomar AO imaging 

Near-infrared adaptive optics imaging of Kepler-14 was 
obtained on the night of 03 July 2010 UT with the 
Palomar Hale 20 0- inch telescope and the PHARO near- 
infrared camera (jHavward et al.ll200lD behin d the Palo- 
mar adaptive optics svstem iTrov 'erall[2000l) . PHARO, 
a 1024 X 1024 HgCdTe infrared array, was utihzed in 25.1 
mas/pixel mode yielding a field of view of 25". Observa- 
tions were performed in both the J (Aq = 1.25 /im) and 
Ks (Aq = 2.145 /Ltm) filters. The data were collected in 
a standard 5-point quincunx dither pattern of 5" steps 
interlaced with an off-source (60" East) sky dither pat- 
tern. The integration time per source was 24.5 seconds 
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Figure 4. MMT/ARIES AO images of Kepler-14 showing a 2.4" X 
2.4" field of view in J and Ks band. The binary nature of Kepler- 
14 is clear from the figure, with the fainter star, B, offset by 0.29 ± 
0.01" to the south-east. 

at J and 12.75 seconds aX Kg- A total of 20 frames 
were acquired at J and Kg for a total on-source inte- 
gration time of 8.5 minutes and 4 minutes, respectively. 
The individual frames were reduced with a set of custom 
IDL routines written for the PHARO camera and were 
combined into a single final image. The adaptive optics 
system guided on the primary target itself; the widths 
of the central cores of the resulting point spread func- 
tions were FWHM = 0.15" at J and FWHM = 0.1" 
at Kg- The final co- added images at J and Kg are shown 
in Figure |3l 

Other than the nearby object, two other sources were 
detected within 10" of the primary target. The second 
closest object is separated from Kepler-14 by 5.5" to the 
northwest and has infrared magnitudes of J = 18.08 ± 
0.04 and Kg = 17.28 ± 0.01. An additional source was 
detected to the southwest at a distance of 6.2" having 
infrared magnitudes of J = 19.00±0.05 and Kg = 18.14± 
0.04. 

The close pair was easily resolved by the adaptive op- 
tics at both J and Kg. The pair is separated by ~ 0.28" 
with a position angle of 142° east of north. The pair 
has magnitude differences of AJ = 0.34 ±0.01 and 
AKs = 0.40 ± 0.01. The brighter infrared (and optical) 
source (component A) is the northwestern star. 

3.3. High-resolution ARIES AO imaging 

High-resolution AO images of Kepler-14 were obtained 
using the ARIES instrument on the 6.5-m MMT. ARIES 
is a near-infrared diffraction-limited imager and spectro- 
graph. On 8 November, 2009 it was operated in the 
f/15 mode, with a 40"x40"field of view and a pixel scale 
of 0.04"/pixel. All images of Kepler-14 had exposure 
times of 10 seconds, with 16 images in J (in a 4-point, 
4"dither pattern) and 19 images taken in Kg (16 in a 
4-point, 4"dither pattern, and 3 images at other offsets). 
The images for each filter were calibrated using standard 
IRAlQ procedures, and combined and sky-subtracted 
using the IRAF task xdimsum. 

In both J and Kg, the binary appearance of Kepler-14 
is clear, with the fainter star, B, offset by 0.29 ± 0.01" 
to the south-east. The separation is comparable to the 
image FWHM (0.5"in J and 0.3"in Kg). The rela- 
tive magnitudes were estimated by PSF fitting, yielding 
A J = 0.398 ± 0.008 and AK = 0.490 ± 0.005. 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



Table 2 

Relative Radial- Velocity Measurements of Kepler-14 



HJD 


Phase 


RV 


i^RV 


BS 


o"BS 


(days) 


(cycles) 


(ms-l) 


(ms"-') 


(m s"-"-) 


(ms"-') 


2455048.454298 


11.394 


-246.1 


25.0 


21.4 


7.4 


2455052.427895 


11.979 


31.5 


17.0 


4.2 


6.5 


2455107.428247 


20.079 


-219.1 


14.2 


-6.8 


4.9 


2455108.417046 


20.225 


-385.2 


15.4 


-4.0 


5.4 


2455109.356436 


20.363 


-285.7 


20.3 


15.0 


10.8 


2455109.415242 


20.372 


-269.0 


20.7 


-1.7 


5.3 


2455111.453147 


20.672 


349.5 


14.2 


-30.5 


6.9 


2455112.462361 


20.821 


393.2 


18.0 


-41.2 


8.7 


2455113.456678 


20.967 


89.2 


19.6 


-5.8 


7.0 


2455114.494316 


21.120 


-260.0 


25.6 


18.4 


6.3 


2455115.492291 


21.267 


-404.0 


30.2 


51.2 


11.0 


2455122.475354 


22.295 


-387.8 


20.3 


-0.3 


11.2 


2455123.417485 


22.434 


-152.6 


19.4 


9.4 


7.3 


2455125.407529 


22.727 


398.5 


19.4 


-29.4 


7.3 
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Figure 5. Upper panel: Radial velocity measurements from the 
Fibre-fed Echelle Spectrograph (FIES) at the 2.5 m Nordic Optical 
Telescope (NOT) at La Palma as a function of orbital phase with 
the best orbital fit overplotted. The velocity of the system has been 
subtracted and the fit assumes a circular orbit, fixing the ephemeris 
to that found by the photometry. Middle panel: Phased residuals 
of the velocities after subtracting the best fit model. The rms 
variation of the residuals is 16.2 ms~^. Bottom panel: Variations 
of the bisector spans from the FIES spectra, with the mean value 
subtracted. 

The delta magnitudes from the Palomar and ARIES 
AO imaging are thus similar, but the difference of 0.06 
magnitude in J and 0.09 magnitude in Kg suggests that 
the accuracy is worse than implied by the formal preci- 
sion. We combined all the Speckle and AO imaging re- 
sults for the assessment of the dilution from the nearby 
companion. 

3.4. High-resolution high SNR spectroscopy 
Spectroscopic observations of Kepler-14 were obtained 
using the Fibre-fed Echelle Spectrograph (FIES) at the 
2.5 m Nordic Optical Telescope ( NOT) at La Palma, 
Spain (iDiupvik Ander"sCTl [20Tol ) as well as HIRES 
(|Vogt et al.lll994D mounted on the Keck I telescope on 
Mauna Kea, Hawaii. We acquired 17 FIES spectra be- 
tween 4 August and 20 October 2009, 3 of which were 
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not used in the analysis because of very low SNR due to 
poor observing conditions. One HIRES template spec- 
trum was also observed on 10 September 2009 and used 
to derive stellar parameters. 

For HIRES, we set the spectrometer slit to 0'.'86, re- 
sulting in a resolving power of A/AA ~ 55,000 with a 
wavelength coverage of ~3800-8000A. We reduced the 
HIRES spectrum followin g a p rocedure based on that 
described bv lButler etHI (|199d). 

For FIES, we used the medium and the high-resolution 
fibers (1'.'3 projected diameter) with resolving powers of 
A/AA « 46,000 and 67,000, respectively, giving a wave- 
length coverage of ~ 3600 — 7400 A. We used the wave- 
length range from approximately ~ 3900 — 5800 A to de- 
termine the radial velocities. The exposure time was ap- 
proximately 60 minutes yielding a SNR from 20 to 65 per 
pixel (SNR of 38 to 120 per resolution element) over the 
wavelength range used. The rather large range in SNR 
is due to the variation in instrumental throughput and 
the stellar flux as a function of wavelength, and the lower 
throughput of the high resolution fiber. 

The FIES spectra were rectified and cross correlated 
using a custom-built pipeline designed to provide pre- 
cise radial velocities for Echelle spectrog raphs. The pro- 
cedure s are described in more detail in iBuchhave et al.l 
((20T0I) . The science exposures were bracketed by two 
thorium argon (ThAr) calibration images taken through 
the same fiber and extracted using the same pipeline 
as the science exposures. The ThAr images were then 
combined to form the basis for the fiducial wavelength 
calibration. Once the spectra had been extracted, a 
cross correlation was performed order by order using the 
strongest exposure as the template. The orders were 
cross correlated using a Fast Fourier Transform (EFT) 
and the cross correlation functions (CCFs) for all the or- 
ders were co-added and fitted with a Gaussian function 
to determine the radial velocity. Uncertainties of the in- 
dividual velocities were estimated by a = RMS{v)/Vn, 
where v is the radial velocity of the individual orders and 
N is the number of orders. 

The light from the fainter, but nearly equal magnitude 
stellar companion (B) dilutes the light of the brighter 
star (A). In Section 15.11 use the photometric centroid 
to determine that it is the brighter star (A) which is 
undergoing transit and thus is the planet hosting star. 
The very small angular separation of 0.29"makes it im- 
possible to separate the two stars on the fiber for the 
spectroscopic observations and it is thus necessary to ac- 
count for the effect of the dilution on the measured radial 
velocities (see Section [Ol and ISTT)) . 

In the observed FIES and HIRES spectra, we did not 
see a composite spectrum in any of the observations. We 
would easily have been able to identify two cross corre- 
lation peaks from a composite spectrum, if the two stars 
did not have nearly equal radial velocities. The combi- 
nation of the small angular separation of the two stars 
and the similar radial velocity makes the probability of 
the two stars being a chance alignment highly unlikely, 
and we therefore conclude that the two stellar compo- 
nents are gravitationally bound in a wide orbit yielding 
an undetectable radial velocity offset between the two 
spectra. 

The radial velocity measurements of the combined light 



of the two components in Kepler- 14 are reported in Ta- 
ble [2] The radial velocities are relative, since they are 
measured relative to the strongest of the observed spec- 
tra adopted as the template. We made a separate esti- 
mate of the systemic velocity (the 7 velocity) by corre- 
lating the observed spectra against the synthetic library 
spectrum best matching the stellar parameters. We took 
the mean of these velocities and subtracted the gravi- 
tational redshift of the Sun (0.636 km s~^), which is not 
included in the calculation of the synthetic library spec- 
tra. We found the mean 7 velocity of Kepler-14b to be 
7 = 6.53 ±0.30 kms"^ 

We fitted a circular orbit to the radial velocities re- 
ported in Table [21 adopting the photometric ephemeris, 
leaving only the orbital semi-amplitude, K, and an arbi- 
trary RV offset as free parameters. A plot of the orbital 
solution is shown in the top panel in Figure [5] with the 
residuals to the fit shown in the middle panel. The orbital 
parameters are listed in Table [31 Allowing the eccen- 
tricity to be a free parameter only reduced the velocity 
residuals by a small amount and yielded an eccentricity 
that was insignificant (e = 0.035 ± 0.020). However, we 
included the eccentricity in the light curve analysis in 
Section TS .41 mainlv to allow for more realistic uncertainty 
estimates of the planetary parameters. 

3.5. Bisector analysis 

We carried out a bisec tor span analysis (jOueloz et al.l 
I200H [Torres et al1l2007[ ) of the FIES spectra to explore 
the possibility that the transit like events are due to an 
eclipsing binary blended with light from a third star. The 
bisector spans are plotted in the bottom panel of Figure 

m 

Since the observed spectrum is a composite of the 
planet hosting star and its fainter companion, which we 
assume is stationary in velocity with respect to the reflex 
motion induced by the planet, we expect to see a slight 
asymmetry in the cross correlation peak as a function 
of phase. The bisector span is defined as the velocity of 
the bisector measured near the bottom of the CCF mi- 
nus the velocity measured near the top, and we would 
thus expect a positive bisector span when the host star 
is moving toward us and a negative bisector span when 
the host star is moving away. In the bottom planet in 
Figure [SJ we see a slight hint of this effect, with the bi- 
sector span being predominantly positive around phase 
0.25 and predominantly negative around phase 0.75. The 
amplitude of the bisector spans is significantly less than 
the radial velocity semi-amplitude and the hint of vari- 
ation is in the expected direction, which supports the 
interpretation that the radial velocity variations are due 
to a planetary companion. 

4. WARM-SPITZER OBSERVATIONS 

Kepler- 14 was o bserved during one transit with Wa rm- 
Spttzer /IRAC (We rner et ed] 120041: iFaiio et al.ll2004[ ) at 
4.5 /im (program ID 60028). The observation occurred 
on UT 2010 August 07 and the visit lasted approximately 
14 h 20 min. The data were gathered in full-frame mode 
(256 X 256 pixels) with an exposure time of 30 s per 
image, which yielded 1700 images. 

The method we used to produce phot ometric t ime se- 
ries from the images is described bv .Desert et al.l ([2009) . 
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It consists of finding the centroid position of the stel- 
lar point spread function (PSF) and performing aper- 
ture photometry using a circular aperture on individual 
exposures. The images used are the Basic Calibrated 
Data (BCD) delivered by the Spitzer archive. These 
files are corrected for dark current, flat-fielding, detec- 
tor non-linearity and converted into flux units. We con- 
vert the pixel intensities to electrons using the informa- 
tion on the detector gain and exposure time provided 
in the FITS headers. This facilitates the evaluation of 
the photometric errors. We extract the UTC-based Ju- 
lian date for each image from the FITS header (keyword 
DATE_OBS) and correct to mid-exposure. We convert 
to TDB-ba sed BJD using t he UTC2BJlH procedure de- 
veloped by lEastman etall (|20iai This program uses 
the JPL Horizons ephemeris to estimate the position of 
Spitzer Space Telescope during the observations. We 
then correct for transient pixels in each individual im- 
age using a 20-point sliding median filter of the pixel 
intensity versus time. To do so, we compare each pixel's 
intensity to the median of the 10 preceding and 10 follow- 
ing exposures at the same pixel position and we replace 
outliers greater than 4 cr with its median value. The 
fraction of pixels we correct is less than 0.06%. The cen- 
troid position of the stellar PSF is determined using a 
DAOPHOT-type Photometry Procedure, GCNTRD, from 
the IDL Astronomy Librar^r"!. We use the APER routine 
to perform aperture photometry with a circular aper- 
ture of variable radius, using radii of 1.5 to 8 pixels, in 
0.5 steps. The propagated uncertainties are derived as a 
function of the aperture radius; we adopt the one which 
provides the smallest errors. We find that the transit 
depths and errors vary only weakly with the aperture ra- 
dius for all the light-curves analyzed in this project. The 
optimal aperture is found to have a radius of 4.0 pix- 
els. We estimate the background by fitting a Gaussian 
to the central region of the histogram of counts from the 
full array. The contribution of the background to the to- 
tal flux from the stars is low for both observations, from 
0.1% to 0.55% depending on the images. Therefore, pho- 
tometric errors are not dominated by fluctuations in the 
background. We used a sliding median filter to select and 
trim outliers in flux and position greater than 5 a. We 
also discarded the first half-hour of observations, which 
are affected by a significant telescope jitter before stabi- 
lization. The final number of photometric measurements 
used is 1570. The raw time series is presented in the top 
panel of Figure[6l We find that the point-to-point scatter 
in the photometry gives a typical signal-to-noise ratio of 
330 per image, which corresponds to 92% of the theoret- 
ical signal-to-noise. Therefore, the noise is dominated by 
Poisson photon noise. 

5. ANALYSIS 

5.1. Centroid shifts 

We use a comparison of the photometric centroid in- 
and out-of-transit to determine which component con- 
tains the transit event. These centroids have been mea- 
sured for quarters 1-6 using two methods: a) a fit of 
the transit model to the whitened row and column cen- 
troid time series, which provides an average offset in row 

littp: //astroutils . astronomy . ohio-state . edu/time/j 
http: //idlastro . gsf c .nasa. gov/homepage .html 



and column for each quarter, and b) centroiding of quar- 
terly average in- and out-of-transit images, where the 
in-transit average is constructed from all in-transit ob- 
servations in a quarter and the out-of-transit average is 
constructed from placing the same number of points on 
each side of each transit event. 

Both methods measured essentially identical centroid 
offsets. These offsets were used to reconstruct the posi- 
tion on the sky of the transiting object, using t h e meth - 
ods described in appendix A of iJenkins et al.l ()2010a[ ). 
The final reconstructed transit source location is then 
the average of the reconstructed transit position over all 
quarters. The distance of this average reconstructed po- 
sition from component A is 0.025 ± 0.024" (1.04 sigma) 
and from component B is 0.251 ± 0.030" (8.33 sigma). 
We conclude that the transiting object is component A. 

5.2. Spectroscopic parameters of the host star 

As noted in Section 13.41 we cannot separate the two 
stellar components on the fiber of the spectrograph and 
we thus observed the light from both stars in the spec- 
tra. As argued in Section we assumed that the two 
stars are physically associated and that they formed to- 
gether at the same time. Since the stars have nearly the 
same temperature due to their position on the H-R dia- 
gram, we concluded that the small magnitude difference 
would result in an insignificant change in the host star 
parameters (see Section [5751 for details). 

We derived stellar atmospheric parameters from both 
the HIRES template spectrum and the high SNR FIES 
spectra used for the orbit determination, which can all 
be used because they are not contaminated by absorption 
from an iodine cell. 

For the HIRES spectrum, we used an analysis 
package known as Spec troscopy Made Easy (SME; 
iValenti fc Piskunovl 119961) ■ along w ith the atomic line 
database of IValent i fc Fisched (pOOSl ). From the HIRES 
spectrum using SME, we found the following parameters: 
effective temperature TeSi, = 6395 ± 60 K, metallicity 
[Fe/H] = +0.12 ± 0.06 dex, projected rotational veloc- 
ity usini = 7.9 ± l.Okms"^, and stellar surface gravity 
\ogg^ = 4.11 ± 0.10 (cgs). 

For the FIE S spectra, we deriv ed stellar parame t ers fo l- 
lowing [Torres et al.l (|2002D and iBuchhave eFall ([2010), 
and in addition we employ a new fitting scheme which 
is currently still under development and being readied 
for publication, allowing us to extract more precise stel- 
lar parameters from the high SNR FIES spectra. We 
mention these values here as a check on the SME val- 
ues, but adopt the SME values because our tools are still 
under development. From the FIES spectra, we found 
effective temperature Tcs^ = 6378 ± 80 K, metallicity 
[Fe/H] — 0.00 ± 0.12 dex, projected rotational velocity 
usini — 10.6 ± l.Okms"^, and stellar surface gravity 
\ogg^ — 4.02 ± 0.21 (cgs). All values are in good agree- 
ment with the values derived from the SME analysis, 
within the uncertainties, except for the value of wsini. 

5.3. Properties of the host star 

Global properties of the star including the mass and 
radius were determined with the help of the stel l ar evo - 
lution models from the series by IGirardi et all (l2000l) . 
Isochrones for a wide range of ages were compared 
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against the effective temperature and metallicity from 
the Keck/HIRES spectra, and the mean steUar den- 
sity, p*, as an indicator of luminosity. If we as- 
sume a circular orbit, then the mean stellar density is 
closely re lated to the norma lized semim ajor axis al Rj, 
(see, e.g., iSeager fc Mallen-Ornelas,200Si : iSozzetti et al.l 
mm . which is one of the parameters solved for in the 
light curve solutions described below in Section and 
is often more accurate than the spectroscopic log 5. In 
practice we used a/ Ri, rather than p^,, and the compari- 
son with the isochrones was coupled with the light curve 
solutions, which were carried out using the Markov Chain 
Monte Carlo technique. Specifically, we derived a distri- 
bution of stellar properties by comparing the isochrones 
with each value in the a/Ri, chains paired with values 
for the temperature and metallicity drawn from Gaus- 
sian distributions centered on the spectroscopically de- 
termined values and their errors. 

The presence of the visual companion detected in our 
high-resolution imaging adds a complication, as the ex- 
tra flux reduces the depth of the transit and affects its 
overall shape in subtle ways, biasing the a/Ri, parame- 
ter. The impact of this extra dilution depends on the 
magnitude difference of the companion in the Kepler 
band (Aiip), which we expect to be close to (but not 
necessarily the same as) the measured magnitude differ- 
ences in other passbands (AF, Ai?, AJ, AJ and /S.Ks). 
We therefore proceeded by iteration, in parallel with the 
light curve solutions. We initially ignored the dilution 
effect on a/R^,, and inferred the absolute magnitude of 
the target in the Kp band from the best-fit isochrone. 
Assuming the companion is physically associated and 
the two stars share the same isochrone, we then deter- 
mined its mass along the isochrone with the condition 
that the magnitude difference in V be exactly equal to 
the measured value. We then read off the AKp value 
directly from the isochrone. We repeated this using each 
of the other magnitude difference measurements (taking 
those in J and Kg from the MMT and Palomar to be 
independent), and we averaged the resulting seven val- 
ues of AKp to obtain 0.45 ± 0.10 mag. With the corre- 
sponding relative flux _Fb / a new light curve solution 
was carried out, leading to an improved a/R-i, distribu- 
tion. This, in turn, was compared once again with the 
isochrones, and led to a slightly revised brightness differ- 
ence of AKp — 0.44 ± 0.10 mag. A further iteration did 
not change this significantly. 

As described in Section 15.21 we have determined the 
host star parameters from the composite spectra of the 
primary star A and the fainter companion B, since it is 
not possible to separate the two stars on the fiber/slit of 
the spectrographs. We estimate that the adopted magni- 
tude difference of AKp = 0.44±0.10 mag does not signif- 
icantly affect the derived spectroscopic stellar properties 
of the host star A. The companion star B is estimated 
to be only 30 K hotter and have a stronger surface grav- 
ity of 0.15 dex compared to the host star. We therefore 
choose to ignore the effect of the dilution on the stellar 
parameters of the host star. 

The resulting properties of the host star are listed in 
Table 131 in which the values correspond to the mode of 
the distributions and the uncertainties reported are the 
68.3% (Icr) confidence limits defined by the 15.8% and 
84.2% percentiles in the cumulative distributions. 



We estimated the distance to Kepler- 14 based on 
isochrones by comparing against th e measured magn i- 
tudes from the Kepler Input Catalog (iBrown et al.l201ll) . 
We fitted the spectral energy distribution with magni- 
tudes for the two stars taken from the Girardi isochrones 
resulting in a distance estimate of 980 pc. For an aver- 
age angular separation of 0.29" the semi-major axis of 
the visual pair is approximately 280 AU, and with mass 
estimates of 1.51 Mq and 1.39 Mq for the two stars, the 
corresponding period is of the order of 2800 years. 

5.4. Light curve analysis 

We modeled the folded transit light curve assuming 
spherical star and planet having radius ratio Rp/Ri,. The 
second star adds its light to the total light curve with the 
observed flux ratio between stars B and A being Fb /Fa- 
The planet was constrained to a circular Keplerian orbit 
parameterized by a period P, a normalized semi-major 
axis distance a/i?* and an inclination to the sky plane i. 

The normalized transit light curve, /(t), was calculated 
to be 



,Ui,U2 



/ 1 + 



Fb 
Fa 



(1) 



where z{t) is the sky-projected separation of the cen- 
ters of the star and planet and A is the fraction of the 
stellar disk bloc ked by the planet, given analytically by 
IMandel &: Agoll ()2002[ ). The limb darkening coefficients 
Ml and U2 parameterize the radial brightness profile, I{r), 
of a star as 

^ = 1 - Ml (1 - vT^) - M2 (1 - VT^) ' (2) 

The continuously defined model, /(i), was numerically 
integrated before being compared with the long cadence 
Kepler light curve. In detail, for each measured time, tj, 
we take nj uniform samples tj^k = tj + kAtj — rint/2, 
separated by Atj = T-mt/nj, over the long cadence in- 
tegration interval of Tint = 29.4 minutes. The flux at 
tj was found by computing the Gaussian quadrature of 
the continuous model fluxes f(tj^k)- In practice, we took 
Uj — 20 for all times. 

We determined the best-fit model to the data by min- 
imizing the goodness-of-fit statistic including Gaus- 
sian penalties to restrict the flux ratio Fb/Fa, ecosw 
and esinw to agree with the observed constraints: 

^2 ^^{T.-Fsf , [Fb/Fa -OMlf , 



0.06P 



(esinw- 0.035)^ (e cos cj - 0.0006)^ 



0.0172 



0.00992 



(3) 



where F^ is the measured flux at time ts and a is the 
expected statistical error in the flux measurements. We 
selected cr = 82 ppm such that the reduced-x^ was unity 
for the best fit solution. 

We determined the posterior probability distribution 
for the fitted parameters by using a Differential Evolu- 
tion Markov chai n Monte Carlo (DE-MCMC) algorithm 
()Ter Braakll2006f) with a Metropolis-Hasting s jump con- 
dition and a jump acceptance probability conditional on 
the likelihood £' cx exp(— x^/2). We computed chains for 
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a parallel population of 90 members through ~2 million 
generations requiring that approximately 25% of jumps 
were accepted on average amongst all population mem- 
bers. The chains were checked for adequate mixing and 
convergence by visual inspection and by observing that 
the number of links was much larger than the autocor- 
relation length (equal to the number of links at which 
the chain autocorrelation drops below one half) for any 
selected parameter. We report the 15.8% and 84.2% val- 
ues of the cumulative distribution for each parameter, 
marginalizing over the remaining parameters. 

5.5. Analysis of the Warm-Spitzer light curves 

We used a transit light curve model multiplied by in- 
strumental decorrelation functions to measure the tran- 
sit parameters and their unc ertainties from the W arm- 
Spitzer data as described in iDesert et al.l ()2011al ) . We 
computed the transit light curves with the IDL tran- 
sit routine DCCULTNL from Mandcl fc Agoll (|2002[ ). This 
function depends on one parameter: the planet-to-star 
radius ratio Rp/Ri,. The orbital semi- major axis to stel- 
lar radius ratio (system scale) a/R^,, the impact param- 
eter 6, and the time of mid transit Tc were fixed to the 
values derived from the Kepler light curve and corrected 
for the dilution (See Table|3l). We assumed that the limb- 
darkening is well- approximated by a non-linear law at 
infrared wavelengths with four coefficients ( Claret 200(J|) 
that we set to their values computed by iSina (2010,) . 

The Spitzer /IRAC photometry is known to be sys- 
tematically affecte d by the so-called pixel-phase effect 
(see e.g.. [Charbon neau et al. 2005; Knutson ct al. 2008). 
This effect is seen as oscillations in the measured fluxes 
with a period of approximately 70 min (period of the 
telescope pointing jitter) and an amplitude of approxi- 
mately 2% peak-to-peak. We decorrelated our signal in 
each channel using a linear function of time for the base- 
line (two parameters) and a quadratic function of the 
PSF position (four parameters) to correct the data for 
each channel. We performed a simultaneous Le venberg- 
Marquardt least-squares fit (jMarkwardtl l2009f ) to the 
data to determine the transit and instrumental model 
parameters (7 in total). The errors on each photomet- 
ric point were assumed to be identical, and were set to 
the RMS of the residuals of the initial best fit obtained. 
To obt ain an estimate o f the correlated and systematic 
errors (|Pont et al.l 120061 ) in our measurements, we used 
the residual permutation boot strap, o r "Prayer Bead", 
method as described in Dese rt et al.l ([^69). In this 
method, the residuals of the initial fit are shifted system- 
atically and sequentially by one frame, and then added to 
the transit light curve model before fitting again. We al- 
lowed asymmetric error bars spanning 34% of the points 
above and below the median of the distributions to de- 
rive t he 1 a unc ertainties for each parameter as described 
in De sert et al .1 (2011b). 

5.6. Interpretation of the Warm-Spitzer observations 

We compute the theoretical dilution factor by extrap- 
olating the i^T^-band measurements to the Spitzer band- 
pass at 4.5 /xm. We estimate that 36% of the photons 
recorded during the observation come from the compan- 
ion star. We conclude that the presence of the con- 
taminating star decreases the effective transit depth of 



Figure 6. VFarm-5p«izer transit light-curve of Kepler- 14 observed 
in the IRAC band-pass at 4.5 /im. Top panel: raw and unbinned 
transit light-curve. The red solid lines correspond to the best fit 
models which include the time and position instrumental decorre- 
lations as well as the model for the planetary transit (see details 
in Sectional. Middle panel: corrected, binned by 25 minutes and 
normalized transit light-curve with the best fit in red. Bottom 
panel: residuals of the data from the best fit. 
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Figure 7. The effect of the dilution of Kepler-14 on the measured 
radial velocities as a function of flux ratio between the companion 
star and the host star. The horizontal dashed line represents the 
observed semi-amplitude of the system and the vertical dashed line 
represents the adopted fiux ratio of the two stars. If the stars have 
the same brightness, the corrected radial velocity is about twice 
the observed and at large magnitude differences, the corrected ve- 
locities approach the observed, as expected. 

Kepler-14 by a factor 0.61. We measure the transit depth 
(limb-darkening removed) of Kepler-14 at 4.5 /im and 
find 17221^38 ppm uncorrected for the dilution. This 
corresponds to i?p/i?* = 0.0415lo;ooi[7. Applying the di- 
lution correction we find Rp/R-i, = 0.05311;q oo2i which is 
consistent with the value derived from the Kepler pho- 
tometry at better than the 2-cr level. Our Spitzer ob- 
servations provide an independent confirmation that the 
transit signal is achromatic, which supports the plane- 
tary nature of Kepler-14b. 

5.7. Dilution effect on the radial velocities 

The measured radial velocities of the host star (A) are 
affected by the light contributed by the companion star 
(B), because the spectrum of B is assumed to be sta- 
tionary in velocity, while the spectrum of A is Doppler 
shifted due to the gravitational pull of the planet. The 
amplitude of the observed radial velocities will thus be 
smaller than if the light from A had not been diluted, 
since the peak of the cross correlation function (CCF) 
from which we derive the radial velocities will be pulled 
toward the stationary CCF of B. 

In order to assess the dilution effect on the radial veloc- 
ities and thus the semi- amplitude of the orbit, we mod- 
eled the effect using the observed spectrum of Kepler-14. 
We shifted the observed spectrum in 50 ms~^ increments 
and co-added the shifted spectrum, representing star A, 
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with the same observed spectrum divided by a constant 
to simulate the stationary companion B at different flux 
ratios of / = Fb/Fa- We analyzed this composite spec- 
trum using the same tools used to extract the radial ve- 
locities for the orbit. 

The relation between the artificially induced radial ve- 
locity shifts, Vim and the resulting "measured" radial 
velocity shifts of the composite spectrum, Vout, is linear 
at given flux ratio, as expected: Vout = afVin, where a/ 
is the slope at a given flux ratio. Vin thus represents the 
true (corrected) radial velocities of the host star {Vcor) 
and Vout represents the observed radial velocities of the 
host star {Vobs)- 

We carried out this analysis at different flux ratios, 
fitting the linear relation between Vin and Vout , thus ob- 
taining the slope a/ at each fiux ratio. We then fitted 
the slopes, a/, themselves as a function of flux ratio with 
a 3rd order polynomial. This enables us to calculate the 
dilution effect for the system at any flux ratio: 

V - - (A) 

y cor — I -(• I -f9 I -fS ^ 

af CQ + cif + C2P + C3f'i 

where / is the flux ratio of the two stars and Ci are the 
polynomial coefflcients of the 3rd order polynomial. The 
dilution effect on the corrected radial velocities as a func- 
tion of flux ratio can be seen in Figure [T] The horizon- 
tal dotted line represents the observed semi-amplitude 
of the system and the vertical dotted line represents the 
adopted flux ratio of the two stars. As a sanity check, 
we see that if the host star and companion have similar 
brightness, the corrected radial velocity is about twice 
the observed radial velocity and at large magnitude dif- 
ferences, the corrected radial velocity approaches the ob- 
served radial velocity, as expected. 

The observed orbital semi-amplitude of Kepler-14 is 
Kobs — 401.7 ± 7.1 ms^^. Since the two stars are 
nearly the same temperature, the dilution changes only 
minutely as a function of wavelength. We thus used the 
magnitude difference of AKp — 0.44 ± 0.10 in all orders, 
and found the corrected semi-amplitude of the orbit to 
be K^or = 682.9 ± 7.1m s-i. 

5.8. Dilution effect on the planetary parameters 

The dilution of the nearly equal magnitude stellar com- 
panion significantly affects the derived planetary param- 
eters of Kepler-14b. The contamination affects the ob- 
served transit light curve depth and therefore the inferred 
radius ratio. In addition, this dilution has a significant 
effect on the light curve profile affecting the inferred geo- 
metric orbital parameters, most notably the normalized 
semi-major axis, a/i?*. If dilution effects are neglected, 
the mean stellar density estimate - which is acutely sen- 
sitive to a/Ri, - used in conjunction with spectroscopic 
stellar constraints will yield significantly inaccurate de- 
rived stellar properties. 

If we assume that the fiux contribution from B is zero 
(i.e., Fb/Fa = and Kobs = 401.7 ± 7.1 ms"!), we find 
that Rp^nocorr — 1 •036lo!Qg4 i?j . Using the derived mag- 
nitude difference AKp = 0.44 ± 0.10, however, we find 
the planetary radius to be Rp = 1.136]'^o!o54 which is 
almost 10% larger. 

As described in Section [5. 71 the orbital semi-amplitude 
is also significantly affected by the dilution. Us- 
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Figure 8. Mass-radius diagram of currently known transiting 
exoplanets. Kepler planets are shown as red squares and planets 
from other surveys are shown as blue circles. The Solar system 
planets are shown as green triangles. The dotted lines are iso- 
density curves (in gcm~^). Kepler-14 is one of the most massive 
transiting exoplanets discovered. 

ing the observed orbital semi-amplitude of Kobs — 
401.7 ± 7.1 ms^^, the uncorrected mass of Kepler-14b 
is Mp^nocorr = 5.14!to ig Mj. After correction for dilu- 
tion, the semi-amplitude increases to Kcorr = 682.9 ± 
7.1 ms~^, which in turn leads to a planetary mass that 
is significantly larger than before (by ^60%): Mp = 
8.40l°:i^Mj. 

The effect of the dilution is much greater on the mass 
than on the radius of the transiting planet. As described 
above, the dilution of the observed transit light curve 
changes not only the depth of the transit, but also the 
light curve profile which in turn affects the inferred stellar 
density estimate. The radius of the planet is thus not 
affected greatly by the dilution, because these two effects 
work against each other. The stellar mass, however, is 
not strongly affected by the dilution and the effect on 
the planetary mass therefore comes almost entirely from 
correction of the orbital semi-amplitude. 

6. DISCUSSION 

We present the discovery of a transiting hot- Jupiter in 
a close visual binary. Had the visual companion not been 
detected, the planetary parameters for Kepler- 14b would 
have been significantly biased. The dilution {AKp = 
0.44 ± 0.10 in the Kepler band) results in a planetary 
mass that, if left uncorrected, is only 60% of the correct 
value, and a planetary radius that is too small by about 
10%. 

The close angular separation of this physically associ- 
ated visual companion makes it essentially undetectable 
spectroscopically: the similar radial velocity as the main 
star means the spectrum is effectively single-lined, and 
the wide orbit (P ~ 2800 yr) implies motion that is slow 
enough that there are no measurable changes in the ve- 
locity of the primary due to this companion. It is only 
with high-resolution imaging that we were able to detect 
it. Many of the over 120 published transiting planets 
and the over 500 published radial velocity planets have 
not been subjected to high-resolution imaging. It is thus 
possible that some of the published exoplanets have in- 
correctly determined planetary parameters, if they have a 
stellar companion like Kepler-14 and the companion has 
not been taken into account. Since many of the published 
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Table 3 

System Parameters for Kepler-14 



Parameter 


Uncorrected 


Corrected 






(Adopted) 



Transit and orbital parameters 

Orbital period P (d)''-'^ 

Midtransit time E (HJD)'' 

Transit duration (days)=''S 0-256lt°-oo25 

Scaled semimajor axis a/R-t'^'^ 7A72tQ''^^\ 

Scaled planet radius Rp/Rt,'^'^ 0-0448^0 0002 

Impact parameter b = acos t//?*^''' 0.531^g Q^g 

Orbital inclination i (deg)''''^ 85.921° *° 

Orbital semi-amplitude K (ms~i)'''= 401.7 ± 7.1 

esiniij^''^ 

ecosoj^''^ 

Center-of-mass velocity 7 (kms"-'-)'^ 

Observed stellar parameters 

Kepler magnitude Kp ' 

Effective temperature T(,ff (K)"^ 

Spectroscopic gravity logg (cgs)"* 

Metallicity [Fe/H]d 

Projected rotation vsini (kms~^)'* 

Derived stellar parameters 

Mass M*(M0)d.e ^■'^O'^tomo 

Radius R^IRq)'^'" 2.358t{j 

Surface gravity log (cgs)^.'' 3.899tJ! g4f 

Luminosity L* (Lq)'^'" S.SOtUr 

Age (Gyr)d.<= 2.0t°oi 

Distance (pc)°'* 

Planetary parameters 

Mass Mp (Afj)='>b><;.d,c 5.141°-^^ 

Radius Rp (Rj, equatorial)'''b,c,d,o 1.036to o84 

Density pp (gcm-^)'^-^-''''^''' S.?!}'^ 



6.7901230 ±0.0000043 
2454971.08737 ± 0.00018 
2591+°-°°''° 

Q 0-1 q+0.578 

O.Z10_Q Qgg 



0.0569 ± 0.0013 

o.ool 

90.0"* 



-,+0.41 
-0.00 
-,+0.0 

682.9 ± 7.1 
0.0350 ± 0.0170 
0.0006 ± 0.0099 

6.53 ±0.30 



12.128 
6395 ± 60 
4.11 ±0.10 
+0.12 ± 0.06 
7.9 ±1.0 



1.512 ±0.043 

2 048+°'-'--'-^ 

n 004+0.028 

•J-»»'*_0.036 

"■^^-0.58 
9 9+0.2 
^•^-0.1 

980 



1.40"* 



-,+0.19 
-0.18 
l.iJD_Q (,54 

7.1 ± 1.1 



Based on the Kepler photometry. 

Based on the dilution by the companion star. 

Based on the FIES radial velocities. 

Based on an SME analysis on the HIRES spectra. 

Based on the Girardi stellar evolution models. 

Based on the Kepler Input Catalog. 

First to fourth contact point. 

The actual orbital period differs fractionally from this value by 2.2 ± 0.1 X 10^ 



result of time dilation for the quoted gamma velocity. 



transiting planets have bright host stars, a campaign to 
gather high-resolution imaging of the host stars could be 
carried out with a modest amount of telescope time. 

In this paper we confirm and characterize the plan- 
etary nature o f Kep ler-14b, also known as KOI-98 in 
IBorucki et all (|2Q11[ ). Kepler-14b has a period of P = 
6.7901230±0.0000043days, amass of Mp = 8.40+!^-^^M. 



and a radius of R„ 



1.136 



+0.073 



p,nocorr 
0.084 ^J- 



— 5. 14]^Q j^g A/j and Rp^nocorr 



0.18 - 

0.054 ^J' yielding a mean 
density of pp = 7.1 ± 1.1 gcm"-^. Not taking the dilution 
into account, the derived mass and radius of the planet 
would be M„ - ^ 1.+0.15 

l.UOO_QQg4 

Kepler- 14b is plotted on a mass-radius diagram in Fig- 
ure [51 which shows all the known transiting exoplanets. 
Kepler- 14b is one of the most massive transiting exoplan- 
ets discovered and is situated in a less dense part of the 
mass-radius diagram together with six other planets of 
similar mass. 
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